Gravitational waves from inspiralling compact binaries: hexagonal template 
placement and its efficiency in detecting physical signals. 
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Matched filtering is used to searcli for gravitational waves emitted by inspiralling compact binaries 
in data from the ground-based interferometers. One of the key aspects of the detection process is 
the design of a template bank that covers the astrophysically pertinent parameter space. In an 
earlier paper, we described a template bank that is based on a square lattice. Although robust, we 
showed that the square placement is over-efficient, with the implication that it is computationally 
more demanding than required. In this paper, we present a template bank based on an hexagonal 
lattice, which size is reduced by 40% with respect to the proposed square placement. We describe 
the practical aspects of the hexagonal template bank implementation, its size, and computational 
cost. We have also performed exhaustive simulations to characterize its efficiency and safeness. We 
show that the bank is adequate to search for a wide variety of binary systems (primordial black 
holes, neutron stars and stellar mass black holes) and in data from both current detectors (initial 
LIGO, Virgo and GEO600) as well as future detectors (advanced LIGO and EGO). Remarkably, 
although our template bank placement uses a metric arising from a particular template family, 
namely stationary phase approximation, we show that it can be used successfully with other template 
families (e.g., Pade resummation and effective one-body approximation). This quality of being 
effective for different template families makes the proposed bank suitable for a search that would 
use several of them in parallel (e.g., in a binary black hole search). The hexagonal template bank 
described in this paper is currently used to search for non-spinning inspiralling compact binaries in 
data from the Laser Interferometer Gravitational- Wave Observatory (LIGO). 

PACS numbers: 02.70.-c, 07.05.Kf, 95.75.-z, 95.85.Sz, 97.80.-d 



I. INTRODUCTION 

Ground-based laser interferometer detectors such as 
Laser Interferometer Gravitational wave Observatory 
(LIGO) [11 or Virgo are expected to detect gravi- 
tational wave (GW)-signals in data that have been, or 
will soon be, collected. The most promising and well- 
understood astrophysical sources of gravitational waves 
are inspiralling compact binaries (ICB) in close orbits [3i] , 
which consist of two compact objects such as primor- 
dial black holes, neutron stars and/or stellar-mass black 
holes. 

Potentiality of a detection verges towards one event 
per year. However, the detection rate strongly depends 
on the ICB coalescence rate [1, @ and the volume of 
universe that detectors can probe. While we cannot influ- 
ence the coalescence rates, we can increase the volume or 
distance at which a signal can be detected, which highly 
depends on (i) the design of the detectors and their sen- 
sitivities, and (ii) on the detection technique that is used. 
Detector sensitivity can be increased most certainly; but 
once data have been recorded, only the deployment of 
an optimal method of detection can ensure the highest 
detection probability, and that is a passport, not only to 
probe the largest volume of Universe possible, but also 
to detect a GW-signal directly for the first time. For- 
tunately enough, altough the two body problem cannot 
be solved exactly in general relativity, post-Newtonian 
(hereafter PN) approximation have been used to obtain 
accurate models of the late-time dynamics of ICB [7|. 
Therefore, we can deploy a matched filtering technique. 



which is an optimal method of detection when the sig- 
nal buried in Gaussian and stationary noise is known 
exactly. The models that we used for detection are also 
called template families. 

The shape of the incoming GW-signals depends on var- 
ious parameters, which are not known a priori (e.g., the 
masses of the two component stars in the case of a search 
for non-spinning binaries). Thus, we have no choice but 
to filter the data through a set of templates, which is 
also called a template bank and must cover the parameter 
space that is astrophysically relevant. Since we cannot 
filter the data through an infinitely large number of tem- 
plates the bank is essentially discrete. Consequently, the 
mismatch between any signal and the nearest template 
in the discrete template bank will cause reduction in the 
signal-to-noise ratio (SNR). Spacing between templates 
must be chosen so as to render acceptable this SNR re- 
duction as well as the computational demand required by 
the cross correlation of the data with the entire discrete 
template bank. As we shall see, the spacing between 
templates is set by specifying a minimal match between 
any signal and the template bank. In practice, template 
families are approximation of the true gravitational wave 
signal, and no true signal will perfectly match any of the 
template families. However, in this paper we shall con- 
sider that template and simulated signal belong to the 
same template family. 

The template bank placement is one of the key aspects 
of the detection process. Nonetheless, its design is not 
unique. There are essentially two types of template bank 
placements. The first one does not assume any knowl- 
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edge on the signal manifold; the second does. The first 
type of placement computes matches between surround- 
ing templates until two templates have a match close to 
the requested minimal match, and computes matches re- 
peatedly over the entire parameter space until it is fully 
populated. Using geometrical considerations, an efficient 
instance of this technique has been developed |p . A sec- 
ond approach, described in various papers [U [HI, fl^ . 
utilizes a metric that is defined on the signal manifold. 
It uses local flatness theorem to place templates at proper 
distances [l3l over the parameter space. We developed a 
template l3ank placement in [T2| that was implemented 
and fully tested within the LIGO algorithm library [l^ . 
This template bank was used in the analysis of data from 
different LIGO science runs [H, d, 113, 111, d ■ We also 
shown that although robust with respect to the require- 
ment (matches should be above the minimal match) , it is 
over-efficient. This result was expected because we used 
a square lattice to place templates over the parameter 
space. 

In this paper, we fully describe and validate a hexag- 
onal template bank placement that is currently used by 
the LIGO scientific collaboration so as to analyze the 
most recent science runs. In Section |lTl we recapitu- 
late some fundamental techniques and notions that are 
needed to describe the bank placement, and previous re- 
sults on the square template bank placement. We also 
provide a framework to validate a template bank. In Sec- 
tion [nTl we describe the algorithm that places templates 
on a hexagonal lattice. In Section HVl we summarize the 
outcome of the simulations performed to test the hexago- 
nal bank. We envisage various parameter spaces that al- 
lows to search for binary neutron star (BNS), primordial 
black hole (PBH), black hole - neutron star (BHNS) and, 
or binary black hole (BBH) signals. We also considered 
design sensitivity curves for the current and advanced 
generation of ground-based detectors. In Section IIVBI 
we show that the proposed hexagonal template bank has 
the required specifications. 

Finally, in addition to the case of a template family 
based on the stationary phase approximation, we also in- 
vestigate in Section IIV CI the possibility to use the same 
hexagonal bank placement with other template families 
including Fade resummation and effective one-body ap- 
proximation. We show that there is no need to construct 
specific template bank for each template family: the pro- 
posed bank can be used for the different families that we 
looked at in this paper. 



II. FORMALISM AND TEMPLATE BANK 
VALIDATION 

Matched filtering and template bank placement use 
formalisms that are summarized in this Section. We also 
review the main results of the square placement, and re- 
capitulate the framework introduced in 12] that allows 
us to validate a template bank. 



A. Signal and Metric 

The matched filtering technique is an optimal method 
to detect a known signal, s(t), that is buried in a station- 
ary and Gaussian noise, n{t) [2^. The method performs 
a correlation of the data x{t) = n{t)+s{t) with a template 
h{t). In this paper, we shall assume that s{t) and h{t) 
are generated with the same model so that a template 
can be an exact copy of the signal. Matched filtering of 
the data x(t) with a template h{t) can be expressed via 
the inner product weighted by the noise power spectral 
sensitivity (FSD), Sh{f), and is given by 

(x,/.) = 2y^ MT) ^'-'^ 

Note that for simplicity, we will ignore the time t within 
the inner product expressions. A template and a signal 
can be normalized according to 



The SNR after fihering by h{t) is 

p=^^ = {x:h). (2.3) 

The simulations that we will perform assume that tem- 
plate and signal are normalized, that is {h,h) — 1, and 
(s, s) = 1. In this paper, we are interested in the fraction 
of the optimal SNR obtained by filtering the signal x{t) 
with a set of template h{t), therefore, we can ignore the 
noise n{t), and {x,h) becomes (s,/i). Strictly speaking, 
(s, h) does not refer to a SNR anymore, but to the am- 
biguity function, which is by definition always less than 
or equal to unity if the two waveforms are normalized. 
In the following, we shall use the notion of match intro- 
duced in [l3|; the match between two templates is the 
inner product between two templates that is maximized 
over the time (using the inverse Fourier transform) and 
the initial orbital phase (using a quadrature matched fil- 
tering). 

The incoming signal has unknown parameters and one 
needs to filter the data through a set of templates, i.e., a 
template bank. The templates are characterized by a set 
of p parameters ^ = 0, 1, ..,p — 1. The templates in 
the bank are copies of the signal corresponding to a set of 
values i?f , i = 0,1, .., Ni,—1, where Nt, is the total number 
of templates. A template bank is optimally designed if 
Nb is minimal and if for any signal there always exists at 
least one template in the bank such that 

minmax(s(/^), /i(i9f )) > MM, (2.4) 

where MM is the minimal match mentioned earlier. Usu- 
ally, in searches for ICB, the value of the minimal match 
is set by the user to 95% or 97%, which corresponds 
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to a decrease in detection rate of 15% and 9%, respec- 
tively. Nevertheless, the minimal match may have a much 
smaller value for the first stage of a hierarchical search 
(e.g., 80%), or for a one-stage search of periodic signals 
(e.g., 70% or lower). 

The distance between two infinitesimally separated 
normalized templates on the signal manifold is given 
by MM 

||/i(i?'^-Fdt?^)-/i(i?'^)|p = IIM'^^II^ 

= g^.d^^'dr, (2.5) 

where /i^ is the partial derivative of the signal h with 
respect to the parameter i?'^. So, the quadratic form 



9^,^ = (/i^t, K) 



(2.6) 



defines the g^^i, metric induced on the signal manifold. 
The metric is used to place templates at equal distance 
in the parameter space. The distance dxi between tem- 
plates in each dimension is given by 



dxi 



MM 



-,z = 0,l,..,j5- 1 



(2.7) 



In practice, using such dxi leaves a fraction of the pa- 
rameter space uncovered, and overlap between templates 
is required (e.g., in the square placement, spacing is ac- 
tually set to dxi/^/2). 

Since we restrict ourself to the case of non-spinning 
waveforms, h{t) depends on 4 parameters only: the two 
component masses, mi and m2 which may vary from 
sub-solar mass to tens of solar mass systems, the initial 
orbital phase (pc, and the time of coalescence tc- We can 
maximized over tc and ipc analytically, therefore the pa- 
rameter space that we need to cover with our template 
bank is a 2-dimensional space only. For conciseness, we 
can represent the GW- waveform with a simplified expres- 
sion given by 



hit) = 



AA-qM 
D 



[T:MJ{t)f'\os[f{t) + ^c], (2.8) 



where f{t) is the (invariant) instantaneous frequency of 
the signal measured by a remote observer, the phase of 
the signal Lp{t) = 2tt J*'^ f{t)dt is defined so that it is 
zero when the binary coalesces at time t = tc, and A 
is a numerical constant representing the amplitude [2]| . 
The asymmetric mass ratio is = mim2/M^, where 
M = mi + m2 is the total mass of the system. There 
exist amplitude corrections up to 2.5PN |8|, the impor- 
tance of which for detection and estimation is shown 
in (2^ . However, in this work, we use restricted post- 
Newtonian models only and limit PN-expansion of the 
phase to 2PN order. Moreover, in the template bank 
placement, namely for the metric computation, we con- 
sider the stationary phase approximation (SP A) 23j . for 
which the metric can be derived analytically [12 1. Nev- 
ertheless, other template families can be used both for 
injection and filtering (see Section flV Ap . 



B. Example : the Square Template Bank 

The placement that we proposed in |12] uses the metric 
based on the SPA model, and the spacing dxi , as defined 
in Eq. (|2.7p . Since the model h explicitly depends on the 
two mass parameters M and 77, then the spacing dxi are 
function of these two quantities as well. However, the 
metric expressed in these two coordinates is not a con- 
stant; it is not a constant either if we were to use the 
component masses, mi and m2. The preference of chirp- 
times, denoted tq and T3 (see appendix |B1 Eqs. IBip as 
coordinates on the signal manifold is indeed more practi- 
cal because these variables are almost Cartesian [23ll23|. 
Although not perfectly constant for PN-order larger than 
IPN, we shall assume that the metric is essentially con- 
stant in the local vicinity of every point on the manifold. 
We could use any combinations of chirptimcs, but using 
the pair tq — T3, there exists analytical inversion with the 
pair M — 1] (see appendix IB2p . 

The parameter space to be covered is defined by the 
minimum and maximum component masses of the sys- 
tems considered (mmin and mmax), and possibly the min- 
imum and maximum total mass (Mmin and Mmax) as 
shown in Fig. [TJ The lower cut-off frequency /l, at 
which the template starts in frequency, sets the length of 
the templates and therefore directly influences the met- 
ric components, the parameter space, and the number 
of templates Nh. In [l2|, we showed how the size of the 
template bank changes with f^. We also investigated the 
loss of match due to the choice of /l . We generally set 
Jl so that the loss of match is of the order of a percent. 

We briefly remind how the proposed square template 
bank works. First, templates are placed along the mi = 
TO2 or 77 = 0.25 line starting from the minimum to the 
maximum mass. Then, additional templates are placed 
so as to cover the remaining part of the parameter space, 
in rows, starting at 77 = 0.25 along lines of constant T3 un- 
til a template lies outside the parameter space. The spac- 
ing between lines is set adequately. Distances between 
templates are based on a square lattice. An example of 
such a placement is shown in Fig. [51 One of the limita- 
tions of the placement is that templates are not placed 
along the eigenvectors of the metric but along the stan- 
dard basis vectors that describe the to,T3 space. This 
approximation make the ellipses slightly more overlap- 
ping than requested and may also create holes when the 
orientation of the ellipses varies significantly (i.e., at high 
mass regime) . The square placement is also over-efficient 
as compared to a hexagonal placement (see Fig. [3]). 



C. Bank Efficiencies 

Independently of the template bank placement, the 
template bank must be validated to check whether it 
fulfills the requirements (e.g., from Eq. 12. 4p . First, we 
perform Monte-Carlo simulations so as to compute the 
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FIG. 1: Example of parameter space and template bank 
placement. The parameter space is defined by the individ- 
ual mass components (from 3 AIq to 30 Mq) and the lower 
cut-off frequency {/l = 40 Hz). The bottom line corre- 
sponds to mi = m2, or equivalently to rj = 0.25. The two 
other boundaries meet where mi — mmin and m2 — mmax. 
The bottom left point of the parameter space corresponds to 
mi — m2 — mmax whereas the top right point corresponds 
to mi — m2 — mmin- The gray/blue circles give the position 
of each template that is needed to cover the entire parameter 
space (black curves). Even though some templates lie outside 
the parameter space boundaries, these are required to fully 
cover the parameter space. 



efficiency vector, £, given by 



£ {Xs,Xh) 



max 

j 



(2.9) 



i=l..N,,j = l..Nt 



where Nb is the number of templates in the bank, Ns the 
number of injections. 

The vectors and correspond to the parameters 
of the simulated signals and the templates, Xs a-nd Xh 
are the models used in the generation of the signal and 
template, respectively. In all the simulations, we set 
1?* — {nil, 1712, Lpc,tc}- Furthermore, we can analyti- 
cally maximize over the unknown orbital phase ipc and, 
therefore, ■d'^ = {7711,7712}. 

The efficiency vector £ and the signal parameter vec- 
tor I?'* are useful to derive several figures of merit. The 
cumulative distribution of £ (Fig. [31 bottom panel) indi- 
cates how quickly matches drop as the minimal match is 
reached. Nevertheless, the cumulative distribution func- 
tion of £ hides the dependency of the matches upon 
masses. Therefore, we also need to look at the distri- 
bution of £ versus total mass M (e.g., Fig. [31 top panel), 
or versus 77, or chirp mass, A4 (see appendix for an exact 
definition). Usually, we look at £]\i only. Indeed, in most 
cases, the dynamical range of rj is small (from 0.1875 to 
0.25 in the BNS case). Finally, we can quantify the ef- 
ficiency of a template bank with a unique value, that is 
the safeness, S, given by 




1.2 ,1.25 , J. 3 
(seconds) 



1.4 




.1 1.15 



1.2 1.25 , 1.3 1.35 

(seconds) 



1.4 
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(2.10) 



FIG. 2: Two instances of template bank placements. In the 
two plots, we focus on a small area of the parameter space pre- 
sented in figure [1] We used a square (top panel) and hexag- 
onal (bottom panel) placement. For convenience, we re-scale 
the metric components so that, goo ~ gn. Each template po- 
sition is represented by a small circle. Around each template 
position, we plot an ellipse that represents iso-match contour 
of MM = 0.95%. Each ellipse contains an inscribed square 
or hexagon which emphasizes how ellipses overlap each other. 
We can see that squares (top) slightly overlap each other. This 
is because templates are layed along ra equal constant line and 
not along the eigen-vector directions, which change over the 
parameter space. In the hexagonal placement, we take care of 
this problem short-coming, and therefore hexagons are per- 
fectly adjacent to each other: the placement is optimal. 



Ideally, we should have a template bank such that S-ji ~ 
MAI. Sti is a generalization of the left hand side of 
Eq. 12.41 on Ng injections. The higher Ng is, the more 
confident we are with the value of the safeness. Ideally, 
the number Ng should be several times the size of the 
template bank that is Ng ^ Nj,, so that statistically we 
have at least one injection per template. The sub-index 
TZ of the safeness is the ratio between Ng and Nj, and 
indicates the relevance of the simulations. The safeness 
provides also a way of characterizing the template bank: 
if S-ji is less than the expected minimal match MM, then 
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the bank is under- efficient. Conversely, a template bank 
can be over-efficient like in Fig. [3l 




0.94 0.96 0.98 1 

l\/latch(%) 



FIG. 3: Efficiencies of the square template bank. For con- 
venience we remind the reader of some results of the square 
template bank provided in [U]. In the simulations, we used 
stationary phase approximant models for both injections and 
templates. Injections consist of binary neutron stars. We used 
4 design sensitivity curves (LIGOI, advanced LIGO, VIRGO 
and GEO), and for each of them we performed 10,000 injec- 
tions. In the top panel, we show all the results together; all 
injections are recovered with a match higher than 95%, as 
requested. In the bottom panel, we decomposed the 4 simu- 
lations and show that all of them behave similarly. Actually, 
we can see that most of the injections are recovered with even 
higher matches (above 97%) showing the over-efBciency of the 
placement. 



III. HEXAGONAL PLACEMENT BASED ON 
THE METRIC 

In the To,T3 basis vectors, both amplitude and orien- 
tation of the eigenvectors change, which may imply a 
laborious placement. In this Section, we describe the 
hexagonal placement that is conceptually different from 



the square placement and takes into account the eigen- 
vectors change throughout the parameter space. 

A. Algorithm 

Although the hexagonal placement algorithm is inde- 
pendent of any genetic or evolutionary algorithms, it can 
be compared to biological process, and we will use this 
analogy to explain the placement. First, let us introduce 
a cell that contain a template position (e.g., tq — T3), the 
metric components defined at this position, and a unique 
identification number that we refer to as an ID. A cell 
covers an area defined by an ellipse with semi-axis equal 
to dxi/2. The goal of a cell is to populate the parame- 
ter space with an offspring of at most 6 cells (hexagonal 
placement). A cell can be characterized by the following 
principles: 

1- Initialization A cell is created at a given position in 

T0 — T3 plane, not necessarily at a physical place (i.e, 
rj can be less than 1 /4) . The initialization requires 
that 

• metric components at (tq — T3) are calculated, 

• a unique ID number is assigned, 

• 6 connectors are created and set to zero. 

Finally, if the cell area intersects with the param- 
eter space, then it has the ability to survive in its 
environment: it is fertile. Conversely, a cell whose 
coverage is entirely outside the parameter space is 
sterile. 

2- Reproduction A fertile cell can reproduce into 6 po- 

sitions that are the corner of a hexagon inscribed 
in the ellipse whose semi-axes are derived from the 
metric components dxi^s. A cell that has repro- 
duced is a mother cell and its offspring is composed 
of 6 daughter cells. Once a daughter cell is initial- 
ized, it cannot reproduce in place of its mother. 
This is taken into account via the connection prin- 
ciple. 

3- Connection Following the reproduction process, a 

mother cell sets the connections with its daugh- 
ter cells by sharing their IDs. Therefore, a mother 
cell knows the IDs of its daughter cells and vice- 
versa. Moreover, when a mother cell reproduces, 
it also sets up the connections between two adja- 
cent daughters so that they both know their IDs. 
These connections prevent cells to reproduce in a 
direction that is already populated. 

4- Sterility A cell becomes sterile (cannot reproduce 

anymore) when both reproduction and connection 
principles have been applied. A cell that is outside 
the parameter space is also sterile (checked during 
the initialization). 
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(a)Initialization (b)Reproduction cell 1 (c) Connections 




(d) Initialization (e)Reproduction cell 2 (f)Reproduction cell 3 



FIG. 4: Hexagonal template bank placement. Using the terminology that is introduced in the text, we can describe the template 
bank placement algorithm as follows. First, in sub-figure (a), a cell/template is arbitrary placed in the parameter space. Its 
coordinates correspond to mi = mmin and m2 = rrimax, and its Id is 1. In sub-figure (b), the cell with ID=1, which is also a 
mother cell, reproduces 6 new cells according to an optimal hexagonal lattice that takes into account its metric components. 
In sub-figure (c) the connection between the offsprings and the mother cell are created (black arrows). Cells that belong to 
the same generation (white ellipses) are also connected if there are adjacent to each other (red arrows). In sub-figure (d), the 
new cells can start to reproduce. However, the reproduction is exclusive: reproduction takes place cell after cell, and the cell 
with lowest ID is chosen to reproduce first. In sub-figure (e), therefore, the cell with ID equals 2 starts to reproduce. Because 
connections already exist with other cells, this cell will reproduce only 3 directions (i.e., 8, 9 and 10). In sub-figure (f), the cell 
with ID equals 3 starts to populate. The 3 new cells (gray ellipses) created by the cell with ID=2 have to wait for the current 
generation (blue ellipse) to fully reproduce. The cells spread until the boundary of the parameter space is reached. 



5- Exclusivity The reproduction process is exclusive: 
only one cell at a time can reproduce. It is exclu- 
sive because a cell cannot start to reproduce while 
another cell is still reproducing. 

The cell population evolves by the reproduction of their 
individuals over as many generations as needed to cover 
the entire parameter space. The first generation is com- 
posed of one cell only. The position of this first cell cor- 
responds to nil = mmin,m2 = mmax- We could start at 
any place in the parameter space. However, local flat- 
ness is an approximation and the author thinks it is bet- 
ter NOT to start at mi ~ rTimax, 1^2 — "^-max where the 
metric evolves quicker (highest mass) . The first cell is ini- 
tiahzed (first principle). Then, the cell reproduces into 
6 directions (second principle). Once the reproduction 
is over, the connectors between the mother cell and its 
daughters are set (third principle) , and finally, the cell be- 
comes sterile (fourth principle). This loop over the first 



cell has created a new generation of 6 cells, and each cell 
will now follow the four principles again. However, the 
new generation of cells will not be able to reproduce in 6 
directions. Indeed, connectors between the first mother 
cell and its daughters have been set, and therefore the 
new cell generation cannot propagate towards the mother 
direction. Furthermore, the 6 new cells have already 2 
other adjacent cells. Therefore, each cell of the second 
generation can reproduce 3 times only. Moreover, some 
of the cells might be outside the parameter space and 
are sterile by definition. Once a new generation has been 
created, the previous generation must contain sterile cells 
only. The algorithm loops over the new generation while 
there exists fertile cells. The first generation is a particu- 
lar case since it contains only one cell. However, the fol- 
lowing generations are not necessarily made of a unique 
cell, and the reproduction warrants a careful procedure: 
the reproduction takes place cell after cell starting from 
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the smallest ID. Moreover, in agreement with the fifth 
principle, the cells of the newest generation wait until 
all the cells of the previous generation have reproduced. 
The reproduction over generations stops once no more 
fertile cells are present within the population. Since the 
parameter space is finite, the reproduction will automati- 
cally stop. Figure[4]illustrates how the first 3 generations 
populate the parameter space. 

Once the reproduction is over, some cells might be 
outside the physical parameter space, or outside the 
mass range requested. An optional final step consists 
in "pushing back" the corresponding cells inside the pa- 
rameter space. First, we can push back the non-physical 
cells only, that is the cells that are below the 77 = 1/4 
line towards the relevant eigen- vector directions onto the 
77 = 1/4 line. Second, there are other cells for which 
mass parameters correspond to physical masses but that 
are outside the parameter space of interest. Nothing pre- 
vents us from pushing these cells back into the parameter 
space as well. This procedure is especially important in 
regions where the masses of the component objects are 
large. Indeed, keeping templates of mass larger than a 
certain value causes problems owing to the fact that the 
search pipeline uses a fixed lower cut-off frequency and 
the waveforms of mass greater than this value cannot be 
generated. In the simulations presented in this paper, 
we move the cells that are below the = 1/4 boundary, 
and keep the cells that are outside the parameter space 
but with 77 > 1/4. Useful equations that characterize the 
boundaries of the parameter space are provided in ap- 
pendix [B] A flow chart of the algorithm is also presented 
in appendix [Cl 

An example of the proposed hexagonal placement is 
shown in Fig.[TJ In this example, the minimum and max- 
imum individual mass component are 3 Mq and 30 M©, 
and the lower cut-off frequency is of 40 Hz. We can see 
that none of the templates are placed below the equal 
mass line whereas some are placed outside the parameter 
space. Figure m gives another placement example. 



B. Size, Gain and Computational Cost 

The ratio of a circle's surface to the area of a square in- 
scribed within this circle is 7ri?^/(2i?^) « 1.57, where R is 
the circle's radius. The ratio of the same circle's surface 
to an inscribed hexagon equals 7ri?^/(3-\/3i?^/2) « 1.21. 
The ratio of the square surface to the hexagon surface is 
therefore about 29%, which means that about 29% less 
templates are needed to cover a given surface when a 
hexagonal lattice is used instead of a square lattice; com- 
putational cost could be reduced by the same amount. 
Tables H] and [ill summarize the sizes of the proposed 
square and hexagonal template bank placements. The 
hexagonal template bank reduces the number of tem- 
plates by about 40% (see Table. [TTT| . This gain is larger 
than the expected 29%, and is related to the fact that we 
take into account the evolution of the metric (orientation 



of cells/ellipses) on the parameter space. 

Computational time required to generate a hexagonal 
bank appears to be smaller than the square bank. In Ta- 
ble IIV| we record the approximate time needed to gen- 
erate each template bank, which is of the order of a few 
seconds even for template banks as large as 100,000 tem- 
plates. It is also interesting to note that most of the 
computational time is spent in the computation of the 
moments (used by the metric space) rather than in the 
placement algorithm. 

The template bank size depends on various parameters 
such as the minimal match and lower cut-off frequency 
that strongly influence the template bank size. Other pa- 
rameters such as the final frequency at which moments 
are computed, or the sampling frequency may also influ- 
ence the bank size. There are also refinements that can 
be made on the placement itself. Two main issues arise 
from our study. First, the hexagonal placement popu- 
lates the entire parameter space. Yet, parameter space 
is not a square but rather a triangular shape. In the 
corner of the parameter space, a hexagonal placement is 
not needed anymore: a single template overlaps with two 
boundary lines. In this case, hexagonal placement can be 
switched to a bisection placement that places templates 
at equal distances from the two boundary lines. A sec- 
ondary issue is that the hexagonal placement is aligned 
along an eigenvector direction. Nothing prevents us to 
place templates along the other eigenvector direction. It 
seems that this choice affects neither the efficiencies nor 
the template bank size significantly. 

TABLE I: Typical square template bank size. We summarize 
the number of templates of typical square template banks. 
We consider several design sensitivity curves such as LIGO, 
VIRGO, . . . (see appendix |X] for analytical expressions and 
lower cut-off frequencies), and 4 typical parameter spaces (see 
Section [TVl for the mass range. 



Bank size 


EGO 


GEO 600 


LIGO-I 


LIGO-A 


Virgo 


BBH 


5582 


1229 


744 


2238 


4413 


BHNS 


94651 


16409 


9964 


35869 


74276 


BNS 


22413 


5317 


3452 


9743 


17764 


PBH 


303168 


62608 


39118 


122995 


242609 



TABLE II: Typical hexagonal template bank size. We sum- 
marize the number of templates of typical hexagonal template 
banks. We consider several design sensitivity curves such as 
LIGO, VIRGO, . . . (see appendix lAl for analytical expressions 
and lower cut-off frequencies) , and 4 typical parameter spaces 
(see Section HVl for the mass range. 



Bank size 


EGO 


GEO 600 


LIGO-I 


LIGO-A 


Virgo 


BBH 


4109 


838 


532 


1712 


3283 


BHNS 


71478 


12382 


7838 


27511 


57557 


BNS 


16036 


3576 


2319 


6969 


12958 


PBH 


205439 


41354 


26732 


84154 


167725 
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TABLE III: Size reduction between the square and hexag- 
onal template banks. We summarize the template bank size 
ratio (in percentage) between the hexagonal and square place- 
ments. The ratios are calculated with the numbers provided 
in Tables. |T] and |lTl For various PSDs and parameter spaces, 
we can see that on average the gain is about 40%. 





EGO 


GEO 600 


LIGO-I 


LIGO-A 


Virgo 


average 


BBH 


36 


47 


40 


31 


34 


37.6 


BHNS 


32 


33 


27 


30 


29 


30.2 


BNS 


40 


49 


49 


40 


37 


43.0 


PBH 


48 


51 


46 


46 


45 


47.2 


average 


39 


45 


40.5 


36.75 


36.25 


39.5 



TABLE IV: Computational cost for different template banks. 
We assume MM = 95%, a Vigo-like PSD with /l = 30 Hz, 
a segment duration of 256s, and a sampling of 4 kHz. Most 
of the computation time is spent in the computation of the 
moments, that depend on the duration of the segment. Using 
short duration vector of a couple of seconds, the computa- 
tional time decreases by about 6 seconds showing that the 
time spent in the placement itself is negligible even for large 
template banks. 







square 


Time(s) 


A'^hoxa 


Time(s) 


0.5 


30 


182136 


25.0 


124652 


9.5 


1 


3 


10187 


7.5 


7251 


6.3 


1 


30 


34095 


9 


24501 


7 


3 


30 


2422 


6.3 


1764 


6.1 



EGO PSD, which we have not used previously, we set 
the lower cut-off frequency /l = 20 Hz. Actually, this 
value can be decreased to about 10 Hz for the BBH case, 
increasing the template bank sizes. 

In all the simulations, we tend to use common pa- 
rameters so as to simplify the interpretation. We use 
a sampling frequency of 4096 Hz over all simulations be- 
cause the last stable orbit /lso — l/(6'^/^27rM) is less 
than the Nyquist frequency of 2048 Hz for most of the 
BBH, BHNS, and BNS signals. The computational time 
is strongly related to the size of the vectors, whose length 
depends on the time duration of the template/signal used 
in our simulations. In order to optimize the computa- 
tional cost, in each search, we extract the longest tem- 
plate duration that we round up to the next power of 2. 
The vector duration is then multiply by 2 for safety. We 
set the minimal match to 95%. We considered 5 types 
of template families that are described later. We can es- 
timate the number of simulations. For instance, using 
Ns injections, with 5 different PSDs, 4 searches (BNS, 
BBH, . . . ), and 5 template families , we have a total of 
A^sX5x4x5 = 100 X Ng injections, which need to be 
filtered through Nj, templates. If we approximate Nf, to 
be 10,000 and Ns to be 10,000 as well, it is clear that 
computational cost is huge. In order to speed up the 
simulations, we chose not to filter signals with all the 
available templates, but only a relevant fraction of them 
around the injected signal; this selection is trivial since 
template and signal are based on the same model. 



IV. SIMULATIONS 

The proposed square and hexagonal template bank 
placements are used to search for various ICB in the 
LIGO and GEO 600 GW-data. They are used to search 
for primordial black holes, binary neutron stars, binary 
black holes and a mix of neutron stars and black holes. 
In the past, the parameter space was split into sub-spaces 
that encompass different astrophysical binary systems 
such as PBH, BNS, BBH, and/or BHNS [H [ll, iB, E, 
[T^. We can filter the data through a unique template 
bank that covers the different types of binaries, however, 
we split the parameter space into the same 4 sub-spaces 
that have been used to validate the square template bank 
placement so that we can compare results together. We 
use the same mass range as in our companion paper, 
that is PBH binaries with component masses in the range 
[O.3-1]M0, BNS [1-3]M0, BBH [3-30]Mq, and BHNS 
with one neutron star with component mass in the range 
[1 — 3]M0 and a black hole with component mass in the 
range [3 — 60] Mq, in which case the template bank must 
cover [4 — 63]Mo. We also use the same PSD by incorpo- 
rating the design sensitivities of current detectors (GEO, 
VIRGO and LIGO-I) and advanced detectors (advanced 
LIGO (or LIGO-A), and EGO). Each of the PSDs has a 
design sensitivity curve, provided in Appendix [X] The 
lower cut-off frequencies are the same as in [l2| and are 
summarized in the appendix as well. In the case of the 



A. Description of the Physical Models 

Theoretical calculations using post-Newtonian approx- 
imation of General Relativity give waveforms as expan- 

1/3 

sions in the orbital velocity v, where v = (27rM/(t)) ' . 
The PN expansions are known up to order in ampli- 
tude and in phase. However, we limit this study to 
restricted post-Newtonian, that is all amplitude correc- 
tions are discarded. Moreover, we expand the flux only 
to 2PN order. The energy function E{v) and the flux 
JF(w) are given by 

E{v) = ENY^^kV^'^^Hy) ^ ^N^^^Fy . (4.1) 

fc j 

We can obtain the phase starting from the kinematic 
equations dt = {dt / dE){dE / dv)dv and d4>/dt = 27r/(t) 
and the change of binding energy = —dE/dt giving a 
phasing formula of the form [25] . 

r-f E'{v) , 

t{v) = Ucf + m / ———dv, 
Jv ^{v) 

^{v) = <^,ef + 2 / v^^dv. (4.2) 

There are different ways in which the above equations 
can be solved. For convenience, we introduce labels so 
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FIG. 5: Cumulative efficiencies of the hexagonal template bank. Both template and signal are based on TaylorF2 model. From 
top left to bottom right (clockwise), injection and template bank cover the PBH binary, ENS, BHNS, and BBH inspiralling 
compact binaries. 



as to refer to different physical template families that are 
used within the gravitational wave community and in our 
simulations. 

TaylorTl If we integrate the equations (14. 2p numeri- 
cally, we obtain the so-called TaylorTl model. If 
instead, we use the P-approximant for the energy 
and flux functions [U, [2y| , then one generates the 
PadcTl model. 

TaylorT2 We can also expand E'{v)/T{v) in a Taylor 
expansion in which case the integrals can be solved 
analytically to obtain the phase in terms of poly- 
nomial expressions as a function of v, which corre- 
sponds to TaylorT2 model [2^. This model is not 
used in this paper but results are very similar to 
the TaylorTS model. 

TaylorTS From TaylorT2, t{v) can be inverted and the 
polynomial expression of v used within the expres- 
sion for 4'{v) to obtain an explicit time-domain 
phasing formula in terms of t. This corresponds 
to the TaylorTS model. 

EOB The non-adiabatic models directly integrate the 
equations of motion (as opposed to using the energy 
balance equation) and there is no implicit conser- 
vation of e nerg y used in the orbital dynamics ap- 
proach [m HI m, nil. The EOB maps the real 



two-body conservative dynamics onto an effective 
one-body problem wherein a test mass moves in an 
effective background metric. 

TaylorF2 The phasing formula is expressed in the 
Fourier domain, and is equivalent to the SPA case 
already mentioned. 



B. SPA Model Results 

First, we validate the hexagonal template bank with a 
model based on the SPA (also labelled TaylorF2), used 
to compute the metric components. We set Xs = Xh = 
TaylorF2, and compute £ and S-jz. We intensively tested 
this bank by setting N^, = 200,000 for each PSD and 
each parameter space considered. Using the template 
bank size from Table |TT1 the ratio 7?, between template 
bank size and number of simulations varies from 1.7 to 
375, which is much larger than unity in agreement with 
discussions that arose in Sec. Ill CI The results are sum- 
marized in Fig. [5] and [51 

In Fig. [5l we notice that the hexagonal bank is ef- 
ficient over the entire range of PBH binary, BNS, and 
BHNS searches. Moreover, the safeness is close to the 
minimal match [S-n € [95%-96%]); by looking at the cu- 
mulative efficiencies, the bank seems to be neither under 
or over-efficient. However, looking more closely at £m 
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FIG. 6: Efficiencies of tfie fiexagonal template bank. Botfi template and signal are based on TaylorF2 model. From top left 
to bottom right (clockwise), injection and template bank cover the PBH binary, BNS, BHNS, and BBH inspiralling compact 
binaries. 



(see Fig. [S]), we can identify a small over-efficient region 
in the BHNS case, where the efficiency is always larger 
than 97% for signals with total mass between [4-20] Mq. 

In the BBH case, the bank is also efficient for the vari- 
ous PSDs with total mass between [6 — 40]A'/q, and simi- 
larly to the BHNS case, it is over-efficient (above 97%) for 
systems with total mass between [8-20] Mq. The bank is 
also under-efficient with matches as low as 93% but for 
very high mass systems above 40 Mq. The match below 
the minimal match are related to the LIGO-I PSD only, 
for which the lower cut-off frequency is 40 Hz. For high 
mass and nearly equal mass systems, the waveforms tend 
to be very short and contain only a few cycles: the metric 
is not a good approximation anymore. It also explains 
the feature seen at high mass, that shows some oscil- 
lations in the matches: a single template matches with 
many different injected signals. One solution to prevent 
matches to drop below the minimal match is to refine 
the grid for high mass range by decreasing the distances 
(i.e, increasing MM) between templates in this part of 
the parameter space. However, the high mass also cor- 
respond to the shortest waveforms which lead to a high 
rate of triggers in real data analysis. Therefore it is ad- 
vised not to over-populate the high mass region. Overall, 
the hexagonal placement has the same behavior as in [12] 
but the bank is not over-efficient anymore in most cases. 



C. Non SPA Model Results 

The square and hexagonal template banks are de- 
signed for TaylorF2 model. Yet, models presented in 
Section HVAI do not differ from each other significantly so 



long as u ^ c, which is the case for PBH, BNS waveforms 
and most of the BHNS and BBH waveforms. Therefore, 
we expect the efficiencies of the template banks to be 
equivalent to the SPA- model results. 

The models used in this Section have the same PN- 
order (i.e., 2PN) as in the TaylorF2 model. The sim- 
ulation parameters are identical except the number of 
simulations that is restricted to Ns = 10, 000 for compu- 
tational reasons. Finally, we tested only the BNS, BHNS 
and BBH searches. The PBH using SPA model being 
sufficient for a detection search. 



1. TaylorTl, TaylorTS, PadeTl 

The TaylorTl, TaylorTS and PadeTl models give very 
similar results that are summarized in the Fig. [71 [8] and 
[51 The safeness is greater than the minimal match for 
the BNS and BHNS searches, for all three waveforms. 
More precisely, S-n ~ 95% for BNS case, and it is slightly 
over-efficient for BHNS case for total mass above 20 M© , 
especially in the case of PadeTl model. In the BBH case, 
the bank is efficient between [6-45] Mq. Then, matches 
drop to 93% for the same reason as in the case of SPA 
discussion. Therefore, we conclude that the proposed 
template bank is also efficient for TaylorTl, TaylorT3 
and PadeTl models. 



2. EOB 

We also investigate the efficiency of the hexagonal 
template bank using EOB templates and signals. The 
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FIG. 7: Hexagonal template bank efficiencies using TaylorTl model. From left to right, results of the BNS, BHNS, and BBH 
injections. 
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FIG. 8: Hexagonal template bank efficiencies using TaylorTS model. From left to right, results of the BNS, BHNS, and BBH 
injections. 



EOB model is intrinsically different from the previous 
models. The results are summarized in Fig. [TOl The 
safeness is slightly under the requested minimal match 
{S-ji = 94.5% w 95%). The template bank is efficient for 
BNS, BHNS and BBH cases. There is no over-efficiency 
noticed in any of the mass range considered. We can also 
notice that the cumulative £m drops quickly and there- 
fore we think that the proposed bank can be used with 



EOB model as well. 



V. DISCUSSION AND CONCLUSIONS 

In this paper, we described a hexagonal template bank 
placement for the search of non-spinning inspiralling 
compact binaries in ground-based interferometers such 
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FIG. 9: Hexagonal template bank efficiencies using PadeTl model. From left ot right, results of the BNS, BHNS, and BBH 
injections. 
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FIG. 10: Hexagonal template bank efficiencies using EOB model. From left ot right, results of the BNS, BHNS, and BBH 
injections. 



as LIGO. The placement is based on a metric computed 
on the signal manifold of a stationary phase approxima- 
tion model. The proposed hexagonal template bank size 
is about 40% smaller than the square template placement 
that was previously used to analyze LIGO science runs 
(i.e., [l^). Yet, the matches between signal and tem- 
plates are above the required minimal match. Therefore, 
the template bank described in this paper is not over- 



efficient: it behaves as required. The main consequence 
is a reduction of 40% of the computational cost required 
to search for inspiralling compact binaries with respect 
to previous searches. 

The proposed template bank is not unique. Several 
parameters can be tuned such as the sampling frequency, 
the final frequency used in the computation of the mo- 
ments, the placement of the template along one eigen- 
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vector or the other, each of which can be investigated in 
more detail. 

The bank was tested with the aid of many simulations 
that use design sensitivity curves for advanced and cur- 
rent detectors, and various inspiralling compact binaries 
with total mass between [O.6-63]M0. We used a model 
based on stationary phase approximation and showed 
that the template bank is efficient for most of the pa- 
rameter space considered. The higher end of the mass 
range was slightly under efficient in the BBH case but 
this is partly related to the shortness of the signal and 
templates considered. 

The proposed template bank can be used for various 
template families, not only the stationary phase approx- 
imation family. In particular, we tested the TaylorTl, 
TaylorT3, PadeTl, and EOB models at 2PN order, that 
have been used for simulated injections in the various 
LIGO science runs. It is interesting to see that the pro- 
posed template bank is efficient for most of the models 
considered in this paper. It is also worth noticing that in 
some cases the template bank is still over-efficient even 
though the bank size is already reduced by 40% (e.g., 
high mass BHNS injections). 

The models that have been investigated in this paper 
are all based on 2PN order, therefore template families 



based on higher PN-order should be investigated. In the 
future, we also plan to consider the case of amplitude 
corrected waveforms. All simulations presented in this 
paper use the same model for both the template and 
signal generation. It would be interested to see how the 
template bank performs when templates are based on one 
model (say, Pade) and the signals are from another (say, 
EOB). 

This hexagonal template bank is currently used within 
the LIGO project to search for non-spinning inspiralling 
compact binaries in the fifth science run. 
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APPENDIX A: DETECTOR'S POWER 
SPECTRAL DENSITIES 



• Finally, the VIRGO PSD is based on data provided 
by J-Y. Vinct and is approximated by 

Shif) = So |(7.87x)-^-« + ^x-' + (1 + x^) I , (A5) 

where 5*0 = 10.2 x lO-'*^ with /o 500 Hz. The 
lower cut-off frequency is — 20 Hz. 



The simulations that we performed use different PSD 
curves that are used to compute the inner products 
(Eq. 12. ip . The different expressions provided uses the 
quantity x = ///o , where / is the frequency and /o is a 
constant. We summarize the different design sensitivity 
curves that have been used in our simulations together 
with the lower cut-off frequency f^: 

• The EGO PSD ^ is given by 
Shif) =S„ {xP' + aixP- + 



a2- 



l + bix + b2x'^ + bsx^ + bix'^ + hr,x^ + b^x^ 



1 -I- ClX 4- C2X^ -I- CzX^ + CiX* 



(Al) 



where Sq = 1.61 x IQ-^^ and /o = 200 Hz. The 
other parameters are : 

= -4.05, p2 = -0.69, ai = 185.62, as = 232.52, 
bi = 31.184, 62 = -64.72, 63 = 52.24, 64 = -42.16, 
65 = 10.17, 56 = 11.53, ci = 13.58, C2 = -36.46, 
C3 = 18.56, and C4 = 27.43. 

The lower cut-off frequency is Jl = 20 Hz. 
• The GEO PSD is given by [s^ 

ShU) - ^0 {lO-4«a;-30 + 34x-i + 



20 



[1 - X- 



X 

2 0.5a;41 



(A2) 



1 + 0.5x2 



where 5*0 = 10"*^ and /o = 150 Hz. The lower 
cut-off frequency is /l = 40 Hz. 



• The LIGO-I PSD [3^ is given by 

Sh{.f) = So{{A.mx)-'>^+ 

0.16a;"'* '^2 ^ 0.52 + 0.32x2} 



(A3) 



where 5*0 = 9 x lO''^'^ ^nd /„ = 150 Hz. The lower 
cut-off frequency is /l = 40 Hz. 

• The advanced LIGO PSD is based on data provided 
in [sO] and given by 



Sh{f) = {a; 



5a; 



111 



l-x^ + 0.5a:'' 
1 + 0.5x2 



(A4) 



where = IQ-''^ and /o = 215 Hz. The lower 
cut-off frequency is Jl = 20 Hz. 



APPENDIX B: PARAMETER SPACE TOOLS 

1. Basic Relations 

Here is a summary of the relationship between individ- 
ual masses mi, 7712, and the two chirptime parameters tq 
and Ta, that are given by 



To 



— — (TTM/i)-'/^ , ^3 = ^ (^M/i) 
2567r/L77 8/l?7 



1 



-2/3 



(Bl) 



where /l is the lower cut-off frequency of the tem- 
plate/signal, M — mi + m2, and rj — 17111712 / . The 
inversion is straightforward; M and rj are given by 



M 



' V-^C-^Y'. (B2) 



327r2/i To ' ' 8/LT3 V 5t3 ; 



It is convenient to introduce the constants Aq and A3 
given by 



. _ 5 A _ 

^0 — ■ ^3 — 



256 (^/l)'/'' 8(7r/L)'/' 



(B3) 



so that Eq. IBll becomes 



ro = ^M-'/', r3 = ^M-2/3. (B4) 
V V 

Finally, the chirp mass, A4, is given by 

M = rj^^^M (B5) 

that allow to to be expressed as a function of chirp mass 
only: 



To 



= AoM-'^^\ 



(B6) 



2. Parameter Space Boundaries relations 

The parameter space is defined by three boundaries 
(see Fig. [T]). On each of these boundaries, we want to 
express T3 as a function of tq. Using lB4l we can eliminate 
M and express T3 as a function of tq and rj: 



T3 



A3 f r]To 
V \Ao 



2/5 



(B7) 
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We can also eliminate rj, and express T3 as a function of 
To and M: 



T3 



A, 
A, 



-ToM. 



(B8) 



The lower boundary corresponds to mi = m2, or rj — 
1/4. Using Eg. IB7[ we can express T3 as a function of tq 
only 



T3 



4.A. 



AAo 



2/5 



(B9) 



J n=i/i 



The second boundary is defined by mi — mmin and 
m2 in [m,„in — m,„ax]- The third boundary is defined 
by mi = minax and m2 in [m,„in - m,„ax]- On those 
two boundaries, we can assume that mi is set to one 
of the extremity of the mass range, denoted m^- Then 



m2 = 
from 



M - me, and t] = (Mme)/(M 



TO = ^ (M)- 



5/3 



Starting 



(BIO) 



we replace 77 by its expression as a function of M and 
me, and obtain after some algebra a cubic equation of 



the form 



' px 







(Bll) 



where x = M^/'^, p = —Ao/{To/me) and q = —mo = 0, 
where mc is either set to mMin or mMax depending on 
which side of the parameter space we are. The solution 
for X is standard and is given by 



27g2 + 




27 


4p3 j ' 


27g2 + 




27 / 



(B12) 



We replace, M 



in Eq. 



to obtain the value of 



T3 on the boundaries when tq is provided. 



APPENDIX C: FLOW CHART OF THE 
HEXAGONAL PLACEMENT ALGORITHM 
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Initialize the ne 

cells 

* 

Connection 



FIG. 11: Flow chart of the hexagonal placement algorithm. 
See the text for detailled description of the initialization, re- 
production, and connection process. 



optional Cieaning : 
Cells outside the parameter space (high 
mass) are push back towards the 
parameter space. 
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